Angular distribution of photoluminescence as a probe of Bose Condensation of 

trapped excitons 
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Recent experiments on two-dimensional exciton systems have shown the excitons collect in shallow 
in-plane traps. We find that Bose condensation in a trap results in a dramatic change of the exciton 
photoluminescence (PL) angular distribution. The long-range coherence of the condensed state gives 
rise to a sharply focussed peak of radiation in the direction normal to the plane. By comparing the 
PL profile with and without Bose Condensation we provide a simple diagnostic for the existence 
of a Bose condensate. The PL peak has strong temperature dependence due to the thermal order 
parameter phase fluctuations across the system. The angular PL distribution can also be used for 
imaging vortices in the trapped condensate. Vortex phase spatial variation leads to destructive 
interference of PL radiation in certain directions, creating nodes in the PL distribution that imprint 
the vortex configuration. 
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The possibility of Bose condensation in two- 
dimensional exciton systems, such as coupled quantum 
wells, has been actively investigated recently [1-6]. The 
interest in these systems is caused by relatively long re- 
combination times that should make it possible to create 
a sufficiently cold system at high density, and thereby 
reach the condensation transition point. Further, their 
two dimensional nature leads to a repulsive interaction, 
disfavouring the formation of biexcitons. It was proposed 
that the condensation of excitons can be facilitated in the 
presence of in-plane traps [7]. Since the temperature of 
the exciton gas depends on the intensity of the laser used 
to pump excitons, it is advantageous if the excitons can 
reach critical density at lower pumping power. By col- 
lecting excitons in a trap, a lower total flux of excitons 
and hence a lower pumping power is required. Small 
traps of micron size arising due to ambient disorder [3] 
(see also [8]), as well as more shallow traps ca. hundred 
microns wide created artificially by stress [9] have been 
considered. 

The density distribution of excitons in such relatively 
shallow traps (of the order of lmeV/fim) changes very 
little through the condensation transition. Thus, in con- 
trast with the cold atom systems [10,11], the direct spa- 
tial imaging of density is not expected to provide dra- 
matic evidence for condensation. The principal mani- 
festations of Bose condensation discussed to date include 
changes in PL spectrum [2,12] and exciton recombination 
rate [1,7]. 

In this work, we demonstrate that PL angular distribu- 
tion can be used as a sensitive probe of Bose condensation 
in a trap. For excitons confined to move in a plane, in 
recombination momentum conservation perpendicular to 
this plane is relaxed. The intensity of radiation at a par- 
ticular angle is therefore controlled by the distribution 
of in plane momenta of the excitons. Bose condensation 



produces large occupation of the state with zero velocity, 
thereby creating a peak in the momentum distribution, 
which gives rise to a sharp peak in the PL angular profile. 
Well below the transition, once the long-range order in 
the condensate extends across the entire exciton cloud in 
trap, the phases of the condensate order parameter be- 
come correlated throughout the cloud. This makes the 
PL radiation from different parts of the cloud fully co- 
herent, and thus focussed in a cone with opening angle 
ca. A ra d/Z, with A rac i the optical wavelength and L the 
cloud size. 

To extract the PL angular profile, we consider dipole 
radiation of excitons, with momentum conservation per- 
pendicular to the plane relaxed [13,14], 



/(k||) ac iV ex (k||)p ph (k||, ko) | (/| ky.p |f) 



(1) 



where ko — 27r/A ra d is photon wavenumber, A r ex (k||) is 
the number of excitons with momentum ky, p p h(k||,fco) 
the photon density of states, and \i) and |/) represent 
the initial and final electronic states. After |k||| is ex- 
pressed in terms of the angle 9 to the plane normal, 
|k||| = fcosin#, the photon density of states and matrix 
element give a factor (1 + cos 2 6) / cos 6. 

The momentum distribution of excitons depends on 
the state of the Bose gas: 
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We model the system as a gas of interacting bosons mov- 
ing in an external trap potential V(r), 



H[i/,) 



d 2 r (3) 



where A is a momentum independent interaction. 
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In the Thomas- Fermi approximation [15], with the sys- 
tem energy being a functional of local density only, one 
can look for the minimum of the energy (3) with the gra- 
dient term omitted. For an harmonic trap V(r) = \ax 2 
containing N particles, this gives a parabolic distribution 
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while V'|r|>fl — 0- The momentum space profile is 

/ sin kR cos fci? N 



iV ex (k) = 8nR 2 N 



\ (kR) 3 (kR) 2 



(4) 



(5) 



(for brevity, from here on we use k = k||). Oscillations 
in the expression (5) result from the sharp edge of the 
Thomas- Fermi distribution (4). 
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FIG. 1. PL angular profile temperature dependence for 
particle number in the trap N = 10 3 (solid lines) and N = 10 2 
(dashed lines). Temperature dependence is calculated for 
the parameters as discussed in the text, with Tbec given by 
Eq.(6). The intensity for N — 10 3 particles is reduced by a 
factor 1/20 to allow comparison. Inset: PL radiation geome- 
try of a 2D exciton cloud. 

The momentum distribution (5) will produce a sharp 
peak in the PL angular profile if R > A ra( j. For the indi- 
rect exciton line A ra( j ~ 800 nm [3] . By either considering 
the capacitance of such a device, or the momentum inde- 
pendent part of the Fourier transform of the interaction, 
we estimate that A = e 2 d/e w 7.0nm 2 eV where d « 5nm 
is the inter-well spacing and e the dielectric constant in 
GaAs. 

A characteristic scale for the trap potential is a w 
10- 12 eVnm~ 2 (i.e. AV = 0.1 meV for Ar = 10 /an, 
[3]). With these values, since the Thomas- Fermi radius 
R ~ 1.5 iV 1 / 4 /im, a peak should be visible even for small 
particle numbers. 

The PL radiation geometry schematic and the peak, 
along with its temperature dependence discussed below, 
are shown in Fig.l. 



The Bose condensation temperature can be e stimated 
from the noninteracting model [3], Tbec = aN/gm, 
with 5 = 4 the exciton degeneracy factor. Alternatively, 
since in a strongly interacting system the density change 
for T < Tbec is small, Tbec can be estimated from the 
degeneracy temperature at the density /9 max of the centre 
of the Thomas-Fermi profile, giving 



Tbec = 



2ttTl 2 
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(for exciton mass m = 0.21 m e , the value m/2Trh 2 = 
0.438 nm _2 eV _1 ). The two figures for T B ec only differ 



by about a factor of 2, because y Xm/2irh 2 w 2. With 

N = 10 4 we obtain Tbec ~ 2 K. (This corresponds to an 
exciton density of ~ 10~ 4 nm~ 2 , comfortably within the 
dilute limit.) 



FIG. 2. The correspondence between vortices and nodes 
in PL angular profile. Left to right: momentum distribution 
(2) for Thomas-Fermi condensate (4) with a single vortex (7) 
located on the x axis at x/R = —0.3, —0.1, 0.1, 0.3. The 
off-center displacement of nodes in k space is 90°-rotated with 
respect to vortex displacement in position space. 

The angular PL profile of the condensate changes in 
an interesting way in the presence of vortices. Here we 
study the momentum distribution for a condensate with 
n vortices, and show that the vortices are 'imaged' by the 
nodes in k-space. Consider the condensate (4) contain- 
ing n vortices. Using the complex variable z = x + iy in 
the plane to write the phase factors due to the vortices 
located at the points Zj, we have 



V>„(r) = "0o(r-)exp 



(7) 



with tjjQ defined by (4) . Vortices give rise to nodes in the 
Fourier transform of ip n (r) and thereby in the momentum 
distribution A r cx (k) (2). There is a one-to-one relation 
between the configuration of the nodes in iV ex (k) and 
the vortex positions. Before discussing the more compli- 
cated case of phase vortices (7), it is instructive to look 
at an example of a condensate function of the form 



ipn(z) = A Y\(z — Zj) exp (—zz/2R 2 



(8) 



describing n 'vortices' with both the amplitude and phase 
varying in space. In this case the configuration of nodes 
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in the momentum distribution is a very simple one. Eval- 
uating the Fourier transform, we represent the polyno- 
mial IX,=i( z — z j) as a differential operator, and obtain 



V> n (k) = 2nA JJ(-2i3fc - Zj) exp (-i? 2 fcfc/2) 
i=i 

n 

= 2 7 rA]^[(i J R 2 fc-z J )exp(-i? 2 fcfc/2) (9) 

4=1 



with the complex variable fc 
has n nodes at the points 



k x + ik y . This expression 



kj — 



-izj/R 2 , j = l, ■ 



(10) 



arranged in exactly the same way as vortices in real space, 
albeit rotated by 90° in a clockwise direction. 

Remarkably, the property that the pattern of the nodes 
in momentum distribution is a replica of the vortex ar- 
rangement, appears to be generic. Although the replica 
is exact only for 'vortices' of the form (8) and Gaussian 
density profiles, all the qualitative features hold also for 
the phase vortices (7). As an illustration, we consider a 
single phase vortex (7) located some distance x away from 
trap center. As can be seen in Fig. 2, in this case there 
is exactly one node in momentum distribution, shifted 
off center by a distance that increases with x. The node 
displacement orientation is at a 90° angle relative to the 
vortex displacement. A similar observation can be made 
for two vortices (Fig. 3a). 

Also, we consider the configuration of nodes in N cx (k) 
corresponding to a vortex array of the form (7). Approx- 
imately, the nodes in this case also form a 90°-rotated 
array, with the distortion being a function of the inter- 
vortex spacing (Fig. 3b, c). 
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FIG. 3. a) PL profile for 2 vortices placed at x/R = ±0.2; 
b), c) PL profile for hexagonal arrays of 7 and 19 vortices. 
The insets show the configuration of vortices in real space. 

Let us now turn to the discussion of finite tempera- 
ture effects. At non-zero temperatures, in the presence 
of density and phase fluctuations, ip = p + 7re lv , we 
construct a perturbation theory in it, tp. As discussed by 
Hohenberg and Fisher [16], who considered the perturba- 
tion expansion for weakly non-ideal Bose gas, this is only 
strictly valid at very low densities, ln(ln(l/na 2 )) ^> 1. 
Their renormalisation group analysis however suggests 
that no other interactions are relevant. 



There are three temperature regimes, a low tempera- 
ture part dominated by phase fluctuations, a free par- 
ticle like regime, and a vortex dominated regime near 
the transition temperature. We consider only the phase 
fluctuations: it will become clear in a moment that the 
temperature dependence of the PL peak is predominantly 
due to the phase fluctuation effects. In this case, equa- 
tion (2) becomes: 



N cx (k)=J2(p(r) P (r')) 1 



/2 p -!-Dr,r' p ik(r-r') 



(11) 



with D rtr ' = (^(tf(r) — <p(r')) j. The phase fluctuations 
are determined by the energy functional 

ffM = / ^(V^(r))Vr (12) 
so that the Green's function G(r,r') — {(p(r)ip(r')} obeys 







V(/)(r)VG(r,r'))=i(r-r') 



For a uniform density profile the solution is [17] 
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(13) 



(14) 



with the thermal length £t = (Xp/^m) 1 ^ 2 /k^T corre- 
sponding to the energy cutoff in the Planck distribution 
of phase fluctuations. 

One can use this result to estimate the fluctuation ef- 
fect for the trapped gas. It is convenient to rewrite the 
logarithm so as to identify a separation dependent and 
an equal point term, 



~(G(r',r') + G(r,r)-2G(r,r')) = 
2irf3p \ \R J V R 



(15) 



We consider in particular the case where r and r' are well 
separated, so G(r, r') is small compared to G(r, r). (The 
typical separation r — r' in Eq.(ll) is of the order of the 
cloud radius R.) In this case 



AWk) 



J Vp(r)e^ G(l%r) e lk ' r d 2 r 



(16) 



Physically, this result means that the phase fluctuations 
are dominated by the equal point contribution, which 
allows one to make a local density approximation in 
Eq.(15) and write G(r, r) = (m/27r/?p(r)) ln(R/£ T ). 

Let us now discuss the phase fluctuations more sys- 
tematically, and verify the local density approximation. 
We shall consider the condensate in an harmonic trap. 
Returning to (13), rescaling r = Rt, and using t as a 
new variable, p(t) — po(l — t 2 ), we write the solution to 
(13) as a sum of angular modes, 
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u(e-e') 



(17) 



The equation for mode I has the form: 
dt 2 



t \l - t 2 )^f + i(l - 3t 2 )% - / 2 (1 - t 2 ). 9 = td(t - t') 



dt 

The substitutions gi(t) = t ±l f(t) and y = t 2 show this is 
the hypergeometric equation [18], the general solution is 

9l {t,t') = l^h+iiU) [t'yh+tity) -*;'/»_,(*>)] (18) 

where /i is an hypergeometric function: 

h±{t) = F(a,b,c;t 2 ) , a + 6 = c=l±Z, a& = ±Z/2 

and t < (t > ) correspond to the smaller (larger) oit,t'. 

As f — » t', the value of G is dominated by the 
large I terms, for which the functions h±(t) tend to 
l/v/2|Z|(l — t 2 ) and so we can consider the terms: 



gi(t + e,t-e) = 



1 



1 - t 2 21 



t 2l {l + 0(e 2 )} + (l-^ 



The second of the terms in brackets gives a divergence as 
e — > 0, of the form 



lim G(i 



•e,t-e) 



-m 



oc 



^/3 Po (l-t 2 )^2/ 

-m /2e 
2n(3p{t) n V~ 



27 



2! 



(19) 



This matches both the nature of the divergence, and the 
pre-factor (—m/2ir(3p) seen in equation (15), and sup- 
ports the local density approximation (16) for a finite 
system, with p replaced by p(v). 

From the asymptotes of the trapped gas Green's func- 
tions (19), cutting the logarithmic divergence with the 
thermal length, the expression for momentum distribu- 
tion due to far separated excitons (16) becomes: 



JVex(kll) 
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(20) 



This expression has the form of replacing the condensate 
density with a "coherent particle density" , 



p(r) -» p(r) 



P( r ) 

PB 
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(21) 



For an harmonic trap, with the numbers discussed above, 

T[K] 
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The temperature dependence of the PL peak is shown in 
Fig.l. Qualitatively, the peak is suppressed at 



T>T*=T BEC /ln{R/£ T ) 



(22) 



i.e. well below the condensation transition. 

In summary, we propose that the angular profile of 
PL can be used as a diagnostic of Bose condensation in 
2D exciton traps. A peak in the emission in a direction 
normal to the 2D plane appears due to long-range phase 
coherence in the condensate, at temperatures for which 
the coherence length exceeds the wavelength of the emit- 
ted radiation. The distribution of PL radiation in the 
peak is sensitive to the presence of phase textures in the 
condensate. In particular, vortices are imaged by the 
nodes in the PL angular profile. 
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